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^ ; Abstract 

■ We introduce a model of a two-core system, based on an equation of the 
| Ginzburg-Landau (GL) type, coupled to another GL equation, which may be 

CO ' linear or nonlinear. One core is active, featuring intrinsic linear gain, while the 

other one is lossy. The difference from previously studied models involving a 
""^ ■ pair of linearly coupled active and passive cores is that the stabilization of the 

system is provided not by a linear diffusion-like term, but rather by a cubic or 
quintic dissipative term in the active core. Physical realizations of the mod- 
els include systems from nonlinear optics (semiconductor waveguides or optical 
cavities), and a double-cigar-shaped Bose-Einstein condensate with a negative 
scattering length, in which the active "cigar" is an atom laser. The replace- 

■ ment of the diffusion term by the nonlinear loss is principally important, as 
diffusion does not occur in these physical media, while nonlinear loss is possi- 
ble. A stability region for solitary pulses is found in the system's parameter 
space by means of direct simulations. One border of the region is also found 
in an analytical form by means of a perturbation theory. Moving pulses are 
studied too. It is concluded that collisions between them are completely elas- 
tic, provided that the relative velocity is not too small. The pulses withstand 
multiple tunneling through potential barriers. Robust quantum-rachet regimes 
of motion of the pulse in a time-periodic asymmetric potential are found as well. 

PACS number: 42. 81. Dp, 03.75.Lm, 42.65.Tg 
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1 Introduction 



Ginzburg-Landau (GL) equations represent a class of universal mathematical 
models which describe pattern formation in various nonlinear media 1 . One 
of the most fundamental types of patterns are solitary pulses (SPs), in loose 
terms often called solitons. In particular, the simplest generic species of the 
GL equations, viz., the cubic complex one, gives rise to a well- known exact SP 
solution However, this solution is unstable (as the equation includes a linear 
gain term, which makes the zero solution unstable, precluding stability of any 
solitary pattern) . Therefore, search for physically relevant models of the GL type 
that give rise to stable pulses has attracted much attention. One possibility is to 
introduce a cubic-quintic GL equation with linear loss and cubic gain, nonlinear 
stability being provided by a quintic loss term. Stable SPs in equations of 
the latter type have been studied in detail |3]- Another model, which finds 
a straightforward physical realization in terms of dual-core nonlinear optical 
fibers, was proposed in Ref. 0]. In this system, one core carries linear gain, 
while the other one is lossy, the corresponding model being based on a system 
of two linearly coupled cubic GL equations (in fact, the one corresponding to 
the lossy core may be a linear equation). 

Detailed investigations have demonstrated that the latter model supports 
stable stationary SPs E] and their bound states in broad parametric 
regions. Stable moving pulses, randomly wandering ones, and breathers (both 
standing and moving) have also been found in this system |5]. The model 
was subsequently generalized to combine it with dispersion-management [5J and 
wavelength-division- multiplexing (i.e., multi-channel) |10j schemes, which opens 
a way for applications to fiber-optic telecommunications. It was also shown that 
a model of the same type can support stable dark solitons 

In all the cases, the stability of SP solutions in this model was provided 
by a linear dispersive-loss term, which accounted for the bandwidth-limited 
character of the gain in the active core (formally, that term is tantamount to 
diffusion). Without this ingredient, the model can generate only unstable SPs. 
However, the natural bandwidth of the optical gain (which is usually provided 
for by the Erbium dopant ^2) is very broad, therefore the limitation of the 
gain bandwidth should be enforced by optical filters specially inserted into the 
system. A related problem is that, in the case of other physical systems which 
may be described by equations of the GL types, such as optical cavities or planar 
waveguides and Bose-Einstein condensates (BECs, see below), the diffusion term 
is not physically possible at all. 

An alternative way to stabilize SPs is to use a nonlinear dissipative term, 
which may be cubic or quintic. In the applications to nonlinear optics, a cubic 
loss term is naturally generated by two-photon absorption, which is a strong 
effect in semiconductor waveguides and semiconductor-doped glasses (see, e.g., 
recent works ^S] and |2]> respectively, and references therein). The investiga- 
tion of this possibility is a subject of the present work. Actually, the system 
with nonlinear loss and without filtering/diffusion is not just another version of 
the above-mentioned dual-core fiber-optic model, but it also directly applies to 
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the description of other physical media, namely, nonlinear optical cavities and 
BECs, as it is explained below. 

In normalized units, coupled equations of the type described above can be 
cast in the form 

Ut = 7 U + ij2U xx + (iai - o- 2 ) \u\ 2 u + iv, (1) 
v t = - (r + ix) v + i"f 2 v xx + iu, (2) 

where, in the case of the optical systems, u and v are amplitudes of electromag- 
netic waves in two cores of the system, the evolutional variable t is either time or 
propagation distance in the optical cavity (depending on the physical formula- 
tion |15p. or the propagation distance in the dual-core optical fiber, and x is the 
transverse coordinate in the cavity (or in a planar waveguide), or the reduced 
time in the application to the fibers. Further, the term with 7 > in Eq. 0) 
accounts for the gain in the active subsystem, the constant of linear coupling 
between the cores is normalized to be 1, V > is the dissipative constant in 
the lossy subsystem, x is a possible frequency- (if t is time) or wavenumber- (if 
t is the propagation distance) mismatch between the cores, and 72, which may 
be assumed positive, is the dispersion/diffraction coefficient (in fact, 72 may be 
given a fixed value by means of obvious rescaling; 72 = 5 will be chosen below, 
as this value is convenient to display numerical results). Finally, 02 and <j\ are, 
respectively, coefficients of the Kerr nonlinearity and nonlinear loss in the active 
core. It is assumed that, in most cases, the field in the lossy core is much weaker 
than in the active one, therefore nonlinear terms in Eq. (J2J) may be neglected 
6 8 , although properties of SPs in the system including nonlinear terms in the 
lossy subsystem are quite similar to those in the system based on Eqs. Q an d 
l|2J) [SJE]. A principal difference of the model based on Eqs. Q and J2J from 
the previously considered ones is that the coefficient in front of the term iu xx 
in Eqs. (JJJ and J2J) is real, while the nonlinear coefficient in Eq. JD is complex; 
previously, exactly the opposite case was considered EllHj • 

In the application to BECs, each equation and (0 may be realized as 
a one-dimensional Gross-Pitaevskii equation for a condensate in a cigar-shaped 
trap, with the linear coupling induced by tunneling between them (see, e.g., Ref. 

The second derivatives in Eqs. and @ are then the kinetic-energy 
terms, with 72 = h/2m, where m is the atomic mass. The linear-loss term in 
Eq. J2J) accounts for evaporation of condensate atoms in the second trap, while 
the linear gain in Eq. Q assumes that the corresponding trap is arranged as an 
atom-wave laser |17| . Further, the coefficient o~\ in Eq. Q is proportional to 
the scattering length of atomic collisions in the BEC gas, and 02 > accounts 
for effective loss due to two-body collisions |18| . 

In fact, in many cases a dominant contribution to the nonlinear loss in BECs 
is due to three-body collisions [T^], hence the corresponding dissipative term is 
quintic, and the accordingly modified equation Q takes the form (again, with 
o 2 > 0) 

Ut — iu + ij 2 u xx + (iai - <T2\u\ 2 ) \u\ 2 u + iv. (3) 
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Below, we will consider two systems, JTJ, © and ©, (J2J with the cubic and 
quintic nonlinear- loss terms, respectively. 

The rest of the paper is organized as follows. In Section 2 we present ana- 
lytical results for SPs obtained by means of a perturbation theory. Numerical 
results are displayed in Section 3, where a parametric region for the existence of 
stable SPs is identified. In Section 4, we consider moving solitons, demonstrating 
that they collide elastically with each other. We also consider a generalization of 
the model including an x-dependent potential, which is relevant for BECs, and 
may be relevant for the case of optical planar waveguides too. In that case, we 
find that moving SPs can coherently tunnel (many times) through a potential 
barrier. The paper is concluded by Section 5. 



2 Analytical results 

Following the lines of Ref. 0], it is possible to develop analysis of SP solutions, 
considering them as weakly perturbed nonlinear Schrodinger (NLS) solitons. 
The linear coupling between the two equations, as well as the gain and loss 
terms, are treated as perturbations, which are ordered so that the conservative 
one, i.e., the linear coupling, is assumed to be a larger perturbation, while the 
nonconservative terms are assumed to be a smaller perturbation. Thus, in the 
lowest-order approximation, that ignores the loss and gain but takes the coupling 
into regard at the first order of the perturbation theory, the SP has the form 

;/ sech [ 4= ) exp (\ 2 t) , v = V(x) exp f Vi ) , (4) 



al \V^J V 2 J V 2 



where it is assumed that a\ is positive, rj is an arbitrary real constant, which 
is an intrinsic parameter of the soliton family, and the real function V(x) is a 
solution to the linear inhomogeneous equation, 



Note that, in this approximation, the linearly coupled equations conserve 
three dynamical invariants: the Hamiltonian and momentum, expressions for 
which will not be used here, and the norm N, which has the physical meaning 
of energy or power in the applications to optics (depending on the particu- 
lar interpretation - it is energy in the fiber, or power in the cavity or planar 
waveguide), or the total number of atoms in the case of BECs, 



N 



u(x)\ 2 + \v(x)\ 2 dx. (6) 



In the next approximation, when the nonconservative perturbations are 
taken into account, we assume that rj may be a slowly varying function of 
time [then, the expression (l/2)r) 2 t for the SP phase in Eq. Q is replaced by 
(1/2) J rj 2 (t)dt\. An evolution equation for 77 (t) can be derived from the balance 



4 



equation for N . In the case of the cubic loss term, which corresponds to Eq. 
{!}, it is 



dN 
~dt 



j\u(x)\ 2 -T\v(x)\ 2 dx~2a 2 I |u(a;)| 4 £fo. (7) 
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In the case of the quintic loss term, corresponding to Eq. 0, |ii(a;)| 4 in the last 
term of Eq. Q is replaced by \u(x)\ . 

By itself, Eq. Q is an exact one. To derive the evolution equation for r](t) 
from it, we substitute the approximation for u, and make use of Eq. 
(which is solved by means of the Fourier transform) . Then, after straightforward 
calculations following the pattern of Rcf. the resultant evolution equation 
can be obtained in an explicit form if x = 0. With the cubic loss term, it takes 
the form 

where C = (tt 2 /6) + C(3) ~ 2.845, and £ is the Riemann's zeta-function. If the 
loss term is quintic, Eq. Q is replaced by 

-Sr' 5 - CT '- 3 < 9 > 

Stationary SP solutions are selected from the continuous soliton family as 
fixed points (FPs) of Eq. |JS} or ©, i.e., as roots of the expression on the right- 
hand side of the equation. The roots can be found in a simple form in the case 
ofEq. ©: 



?7 F p = (16er 2 )" 



15a 2 -f ± J 15of (15(7^2 _ 16Ca 2 r) 



(10) 



An elementary consideration demonstrates that, in either model, there may 
exist two physical FPs or none, depending on values of the parameters: in the 
case of Eq. JSJ, physical solutions exist if 

2cr^ 7 3 > iCa 2 T, (11) 

and in the case of Eq. 0, this condition is replaced by [as it immediately 
follows from Eq. (fTTTfl ] 

15o-i7 2 > 16Ccr 2 r ; (12) 

note that 72 does not appear in Eqs. i|ll|) an d i|12|) . Further, it follows from 
Eqs. © and 10 that the FP corresponding to a smaller value of rj [for instance, 
the one with the lower sign in Eq. I|10fl ] is unstable, while the FP corresponding 
to larger r\ is stable (within the framework of the present approximation). In 
the next section, these analytical predictions will be compared with results of 
direct numerical simulations. 

Lastly, we notice that the stability of the zero solution is a necessary condi- 
tion for the full stability of any solitary pattern. In turn, simple necessary (but, 
generally speaking, not sufficient) conditions for the zero-solution stability take 
a simple form in the case x = 0|2]: 7<r<l/ 7 (which implies that 7 must 
be smaller than 1). 
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3 Numerical results 




Figure 1: Examples of stable stationary solitary-pulse solutions to (a) Eqs. Q, 
©, and (b) Eqs. ©, ©, found for (a) 7 = 0A,a t =l,a 2 = 0.1, 72 = 5,T = 
1,X = 0, and (b) 7 = 0.5, T = l,j 2 = 5,cti = l,cr 2 = 0.05, x = 0. The 
continuous and dashed curves show and |«(a?)|, respectively. 

In the numerical investigation, SP solutions were looked for as solutions to 
the stationary version of Eqs. i ffl l. PJl or iffil. PJl: then, stability of the obtained 
solutions was tested in simulations of the full time-dependent equations. Typical 
examples of thus found stable SPs in both models are displayed in Fig. 1. 

Results of systematic numerical simulations of Eqs. Q), J5J and @, J5J) 
are summarized in Fig. 2, which displays generic examples of stability regions 
for SPs in the parameter plane ((71,7). These two parameters are chosen for 
the variation, while the others are fixed, as they can be readily varied in the 
experimental studies of BECs: the gain 7 is controlled by intrinsic parameters 
of the atom laser, and the nonlinearity <j\ can be changed by means of the 
Feshbach resonance |20|. In the applications to optics, the gain parameter is 
also easy to vary. 

The analytical predictions 1)11(1 and 1(12(1 for the existence of SPs correlate 
with lower borders of the numerically found stability regions in Figs. 2(a) and 
2(b). As is seen, the perturbative result is accurate enough for relatively large 
values of o\ and small values of 7, when the gain and loss may indeed be 
regarded as small perturbations; with the decrease of o\ and increase of 7, 
the perturbations are no longer small, which explains discrepancy between the 
analytically predicted and numerically found lower borders in Figs. 2(a) and 
2(b) in this case. 

Below the existence border, the perturbation theory predicts that no station- 
ary SP exists; in accord with this, the numerical simulations show that, beneath 
the lower border, any initial pulse decays to zero. On the other hand, the upper 
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Figure 2: A generic example of the stability region for solitary pulses in the 
model (pj, (J2J (a) and (QJ), © (b), as found in a numerical form. The fixed 
parameters are a% — 0.1, 72 = 5, T = l,x = in (a), and 02 = 0.05,72 = 
5,r=l,x = 0in (b). The curve denoted by crosses is a border above which 
some perturbations with high wavenumbers grow and SPs are destroyed. The 
curve denoted by diamonds is a border below which SPs decays to zero. The 
dashed curve in each panel is the existence border for the pulses predicted in 
the analytical form by (a) Eq. (11) or (b) Eq. (12). 



border of the stability regions in Fig. 2 actually bounds not the existence, but 
rather stability of the SPs. Above the upper border, perturbations with high 
wavenumbers grow and destroy the SPs. 

The unperturbed NLS equation supports (bright) solitons only in the case 

72^1 > 0, (13) 

see Eqs. Q and ©. In terms of nonlinear optics, this condition corresponds to 
a combination of spatial diffraction or anomalous temporal dispersion and self- 
focusing nonlinearity, or normal temporal dispersion and self-defocusing nonlin- 
earity; in terms if BECs, it corresponds to the case of negative scattering length 

An issue of considerable interest is whether the addition of dissipative terms 
makes it possible to relax the condition 1)13(1. Note that the exact SP solution 
to the cubic complex GL equation exists irrespective of this condition [2], but 
that solution is always unstable. In Ref. [H], stable SPs were found, in the 
two-core model including the filtering term, for both signs of the product 72Ci- 
However, our result for the present model, in which the filtering term is replaced 
by the cubic or quintic loss, is that the condition (|13|) remains necessary for the 
existence of SPs. In fact, fixing 72 to be a positive constant, it is interesting to 
find the smallest value of the nonlinearity coefficient o~\ up to which the stable 
pulse persists. The result [obtained for Eqs. Q and (J2J, with the cubic loss 
term] is shown in Fig. 3, which displays the smallest SP-supporting value of 
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Figure 3: The minimum value of the nonlinear coefficient a\, up to which stable 
solitary pulses were found in Eqs. QJ, 0- varying the gain 7, as a function of 
the linear-loss coefficient T. The other coefficients were fixed: <72 = 0.1, 72 = 5, 
and \ = 0. 

o\ versus the linear-loss factor T, as found by varying the linear gain 7, while 
(72, 72 and x — were fixed. 

Similar analysis was performed for models differing from the ones considered 
above by adding to Eq. J5J the conservative cubic term itri|w| 2 ti, i.e., essentially 
the same one as in Eq. or ©. The result is that the stability regions for 
SPs are similar to, although somewhat smaller than, those shown in Fig. 2. 

4 Moving pulses and coherent tunneling 

One of principal distinctions between the models with the nonlinear loss and 
ones with the linear diffusion (filtering) is that the models considered in the 
present work share the Galilean invariance with the NLS equation. Due to this 
reason, solutions for SPs moving at an arbitrary velocity can be generated by 
the Galilean transform from any quiescent SP. The moving pulse solution with 
an arbitrary wavenumber k can be expressed as 

u k = u (x - ct, t)e ikx ~ iut , v k = v (x - ct, t)e lkx ~ wt , (14) 

where uo(x,t) and vo(x,t) represent a stationary-SP solution to Eqs. JQ) and 
J2J), and 

c = 272/c,w = 7 2 £: 2 . (15) 

In the application to BECs, where 72 = h/2m, Eqs. (|15|) are tantamount to the 
usual relations for the momentum and kinetic energy of a quantum particle, 

hk = mc = P,fruj = (hk) 2 /(2m) = E kin . (16) 

The availability of the moving pulses suggests numerical experiments aimed 
at simulation of collisions between them. A typical example of the collisions 
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Figure 4: A typical example of multiple collisions between two moving solitons 
in Eqs. Q and J2J with periodic boundary conditions. The parameters are 
o"2 = 0.1,72 = 5,r = l,x = and the spatial period is L = 60. The initial 
wavenumbers of the pulses [see Eqs. Q and are fc = ±0.5. 



is displayed in Fig. 4, that presents results of simulations of Eqs. and J5J 
in a domain with periodic boundary conditions, which gives rise to periodic 
recurrence of the collision. It is obvious from this figure, and can be observed 
as a generic situation, unless the relative velocity of the colliding solitons is 
too small, that the collisions are, practically, completely elastic (if the relative 
velocity is too small, both SPs decay to zero state after the collision) . The same 
is true for the model with the quintic loss. 

In the case of BECs, an external potential, such as a harmonic magnetic 
trap, is usually an important ingredient of the model (in the optical model 
describing the spatial evolution of the fields in coupled planar waveguides, a 
similar term may describe a spatially modulated profile of the refractive index 
in the waveguides). Equations (1) and (2) in the case of BECs confined by the 
external potential U(x) are modified as 

ut = 7it + i(h/2m)u xx + (iai — 02) \u\ 2 u + iv — ifi 1 U{x)u, (17) 
vt = — (r + ix) v + i(h/2m)v xx + iu — ih~ 1 U(x)v. (18) 

We have checked that, in the model introduced above, the addition of the har- 
monic potential U(x) = (K/2) (x — L/2) 2 , where x = L/2 is the central point 
of the magnetic trap, gives rise to very persistent periodic oscillations of the 
SP. Moreover, it can be verified that, in this case, the SP behaves as a perfect 
quasi-particle, obeying the corresponding equation of motion 



d 2 X p dU 



ni- 



di 2 dx 



= -K(X p -L/2), (19) 
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Figure 5: Time evolution of the peak position of the oscillating pulse in the 
harmonic potential U(x) = {K/2)(x - L/2) 2 with K = 0.01 for 7 = 0.4, <j 2 = 
0.1, h= 1, m= 0.1, T = l,x = 0. 



where X p is the coordinate of the pulse's peak (provided that potential's strength 
K is not too large). 

The simplest way to derive Eq. (|19f) is to assume the solution in the same 
boosted form as given by the Galilean transformation l|14f> . (|15|) . but assuming 
that the velocity c may be a slowly varying function of time, i.e., 



// n \ x — J c(t)dt,t j exp 
vq ( x — / c(t)dt, t I exp 



imc(t) 
imc(t) 



(20) 



Then, the approximation l|20|l should be substituted into the balance equation 
for the net field momentum P, which in the presence of the external potential 
has the form 



dP d f +oc 

— = — (-ih)(u*u x + v*v x )dx/ 
dt at /_„ 







+ OC 



(\u\ 2 + \v\ 2 )dx 



U(x)^-(\u\ 2 + \v\ 2 )dx/ j (\u\ 2 + \v\ 2 )dx 



— OO 



(21) 



[note that this definition of the field momentum agrees with the relation mc = P 
in Eq. (0 and f^{\u\ 2 + \v\ 2 )dx is constant in time for Eq. (J2DJ]. Finally, 
identifying X p = J c(t)dt and assuming that the potential U(x) varies on a scale 
which is essentially larger than the internal size of the pulse, one can easily derive 
Eq. I|19l) from Eq. I|21l) by means of the integration by parts. 

Figure 5 displays the time dependence of X p in the harmonic potential 
U(x) = 0.005(:e - L/2) 2 for m = 0.1. The initial peak position is X p (0) = 
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Figure 6: An example of multiple coherent tunneling of a traveling solitary pulse 
across a narrow potential barrier (19) in the model (|17fl . i|18|) . The parameters 
are 7 = 0.4, 02 = 0.1, h = 1, m = 0.1, T =1, x = 0, and the initial 
wavenumber pushing the pulse is k = 0.4. 



L/4 = 15, and the initial wavenumber is k = [see Eqs. (|14|) ]. The time 
evolution of the peak position almost exactly follows the corresponding solution 
X p {t) — 30 — 15 cos(VO.lt) of Eq. (jl9|l . If two pulses are originally placed in 
the trap, they perform persistent oscillations, periodically passing through each 
other in the elastic fashion, like in Fig. 4. 

Another physically relevant possibility (in the application to BECs) is to 
consider a system with a narrow potential wall which separates a broad but 
finite trap into two compartments. In other words, it is a double-well configu- 
ration, which is frequently considered in the context of one-dimensional BECs, 
but usually for the case of positive scattering length (i.e., the self-repulsive non- 
linearity), see, e.g., Ref. .21 . To simulate this situation, the potential U{x) was 
taken (for instance) as 

C 2, \x- 30| > 15 
U(x) = l 1, |x-30|<0.5 . (22) 

[ 0, 0.5 < |x-30| < 15 

In this case, the initial SP was given a velocity corresponding to the kinetic 
energy = 0.8 [see Eq. <|16|l ]. which is lower than the central potential 

barrier in the expression (|22|l . 

The result is displayed in Fig. 6, which demonstrates that the moving pulse 
can tunnel across the central potential wall several times. Note that, unlike 
the case displayed above for the case of the smooth harmonic potential, in 
the present situation the pulse's motion does not obey the classical Newton's 
equation of motion, which may be explained by the fact that the potential 
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Figure 7: An example of the quantum-ratchet effect for traveling pulses in the 
time-periodic sawtooth potential (20) for 7 = 0.4, 02= 0.1, % = 1, m = 
0.1, T = l,x = 0. The initial velocity is 4 in (a) and —4 in (b). The temporal 
dependences of the peak position X p are displayed in both panels. 



function (1221) varies steeply in space, breaking the applicability condition for 
Eq. (|19fl . Further, we notice that, in the course of the multiple tunnelings, the 
norm of the pulse remains practically constant, while its kinetic energy gradually 
decreases. Eventually, the tunneling ceases, and the SP finds itself trapped in 
one compartment. 

If the potential function is asymmetric, spatially periodic, and, simultane- 
ously, time-periodic, the quantum ratchet effect may be observed [221 ? 123 To 
demonstrate this possibility, we adopt a time-periodic sawtooth potential: 



U{x,t) = [l + 0.8cos(27rf)] x 



7(10 -x), 9.5<x<10, 

7(25 -x), 24.5<x<25, 

7(40 -x), 39.5<x<40, (23) 

7(55 -x), 54.5 < x < 55, 

0, otherwise , 



so that the peak amplitude of the potential oscillates as 3.5 [1 + 0.8cos(27rf)]. 
We also assume periodic boundary conditions, the spatial period being 15. 

Figure 7(a) displays the time evolution of the peak position for the SP with 
the initial velocity cq = 4 [see Eq. l|15fl ]. As is seen from the figure, the pulse 
steadily moves in the right direction. Note that the initial kinetic energy of 
the pulse, -Ekin = 0.8, is lower than the average height, ([/ max ) = 3.5, of the 
potential peak, hence progressive motion is possible due to the tunnel effect. 
The kinetic energy decreases as a result of the multiple tunnelings as in Fig. 6; 
however, energy supply is possible in the time-periodic potential. Thus, the 
nearly steady propagation is possible, as is seen in Fig. 7(a). On the other 
hand, if the initial velocity is cq = —4, simulations demonstrate that, while the 
kinetic energy again decreases in time, the energy supply is not efficient enough 
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in this case to compensate the loss, and the pulse gets trapped in a potential 
well after a transient, see Fig. 7(b). Thus, the time-periodic sawtooth potential 
admits only the unidirectional drift to the right, which conforms to the definition 
of quantum rachets |23| . 

5 Conclusion 

In this work, we have proposed a model of a two-core system, based on an 
equation of the Ginzburg-Landau (GL) type, coupled to another GL equation, 
which may be linear or nonlinear. One core is active, being equipped with linear 
gain, while the other one is lossy. The difference from previously considered 
models is that the overall stabilization of the system is provided not by the linear 
filtering (diffusion) term, but rather by cubic or quintic dissipation in the active 
core. Physical realizations of the model include several systems from nonlinear 
optics (semiconductor waveguides or optical cavities), and a double-cigar-shapcd 
BEC, in which one "cigar" is actually an atom laser. The replacement of the 
diffusion term by the nonlinear loss is principally important, as diffusion is 
not possible in these physical systems, while the nonlinear loss may occur. A 
stability region for solitary pulses was identified in the relevant parameter plane 
by means of numerical simulations. One border of the region can be predicted in 
an analytical form by the perturbation theory. Moving pulses were considered 
too, with the conclusion that collisions between them are completely elastic 
(unless the relative velocity is too small), and they withstand multiple tunneling 
through potential barriers without losing their coherence. The existence of the 
robust quantum-rachet regime of motion for the pulses was demonstrated as 
well. 
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